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ABSTRACT 


An  exact  Markov  approach  is  used  to  caicuiate  the 
attrition  suffered  by  a  number  of  identical  targets  subjected 
to  sequential  missile  attacks.  Various  missile-allocation 
situations  are  examined.  An  APL  program  is  presented. 
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A  MARKOV  APPROACH  TO  LARGE  MISSILE  ATTACKS 
Maurice  H,  Mi zrahi 


INTRODUCTION 

This  paper  exactly  calculates  the  attrition  suffered  by  a  number  of 
targets  in  a  formation  subjected  to  sequential  missile  attacks*.  The  targets 
could  be  a  raid  of  enemy  bomber  aircraft  attempting  to  penetrate  several 
layers  of  missile  defenses,  such  as  mi ss i 1 e-f i r i ng  fighter  aircraft  or  surface- 
to-ai r  missi les. 

Specifically,  we  will  calculate  the  distribution  of  the  number  of  sur¬ 
vivors  of  a  set  number  of  attacks  with  given  parameters,  for  various  missile- 
allocation  situations,  and  the  expected  number  of  missiles  fired.  The  emphasis 
will  be  on  eliminating  the  complexity  arising  from  a  large  number  of  missiles 

attacking  simultaneously.  APL  programs  for  these  calculations,  written  by  Paul 
E.  Klebe  and  the  author,  are  presented  in  Appendix  A. 

ASSUMPTIONS  AND  METHODOLOGY 

The  targets  are  identical,  and  the  missiles  in  a  given  attack  are 
identical  and  fired  simultaneously.  Let 

n  be  the  initial  number  of  targets 

M.  be  the  maximum  number  of  missiles  fired  in  the  ith  attack 

p.  be  the  probability  of  kill  of  1  missile  against  1  target 
1  in  the  idh  attack 

k  be  the  total  number  of  attacks. 

Three  cases  will  be  examined,  differing  in  the  manner  in  which  the  missiles 
allocate  themselves  and  the  number  of  missiles  fired.  The  first  2  cases  reflect 
random  targeting:  the  missiles  allocate  themselves  randomly  among  all  the 
available  targets  (Case  l)  or  among  a  specific  number  b.  of  targets  --  say  the 
closest  b.  targets  (Case  II).  Under  those  assumptions,  if  no  target  is  present 
no  missile  is  launched,  but  if  at  least  1  target  is  present  all  M.  missiles  are 
launched.  This  corresponds  to  the  situation  where  the  defense  sees  only  a  blob 


of  targets  and  cannot  tell  them  apart,  or  where  the  missiles  cannot  be  pre¬ 
cisely  guided*  The  last  case  (Case  111)  reflects  perfect  targeting:  the 
targets  can  be  distinguished  individually  and  only  1  missile  is  allocated 
per  target*  up  to  a  total  of  M.  in  the  i_th  attack.  Thus,  no  target  can  be 
hit  by  more  than  1  missile  in  Case  1  1 1 u 

It  is  assumed  that  the  targets  do  not  shoot,  back,  or  that,  if  they  do, 
an  adjustment  is  made  to  the  M.  of  the  next  attack.  A  target  is  either  killed 

fl 

(it  "disappears")  or  it  is  intact:  no  intermediate  states  are  considered. 


A  simple  Markov  chain  will  be  used  in  every  case.  The  state  of  the  system 
after  the  itf[  attack  is  therefore  described  by  an  (n-H)  -  dimensional  vector, 


S. 

i 


<po  p, 


e  ) 

n' 


(1) 


where  is  the  probability  that  r  targets  are  present  (have  survived).  The 

initial  state  is 


SQ  =  (0  0  ...  Q  1)  .  (2) 

The  transition  matrix  T1  corresponding  to  the  itji  attack  depends  on  p.  and 
M.  (and  also  on  b.  in  Case  II).  T1  is  an  (n+l)  by  (n-s-l),  triangular, 

(M.  +  1)  -  diagonal  matrix.  After  the  kth  attack,  the  system  is  in  the  state 

Virk  -  so T' T*  ••• T“  •  <3> 


The  problem,  therefore,  reduces  to  calculating  the  transition  matrices 
corresponding  to  each  of  the  3  cases.  Once  that  is  done  and  the  final  state 
is  obtained,  one  can  extract  the  relevant  statistics.  Let  S^O)  (for 
i  =  0,1,... n)  denote  the  ij^h  element  of  S^,  Then: 
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The  expected  number  of  targets  surviving  the  k  attacks  is: 


j  sk(j) 


The  standard  deviation  is: 


j2  sk(j)  - 


J  sk(j) 


®  The  expected  number  of  missiles  fired  is: 


Ek  ■>/<')  - 


where  E(i)  is  the  expected  number  of  missiles  fired  on  the  i_th  attack.  In 
Cases  I  and  I  I ,  we  have: 


E ( i )  =  Mj  [l  -  S._,  (0)]  . 


The  above  reads:  M.  missiles  (the  maximum  for  the  ith^  attack)  multiplied  by 
the  probability  that  there  is  at  least  one  target  present  after  the  (i  -  i)  th 
attack. 


i n  Case  III,  we  have: 


E(i)  = 


J  s;_| (j)  for  H.  >  n 


J  s j (J) 


+  M.  +  ...  +  Sj_|(n)J  for  M.  <  n. 
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The  above  reads:  if  the  number  of  targets  present  after  the  (i  -  l)th  attack 
is  less  than  M. ,  then  only  as  many  missiles  as  there  are  targets  are  fired, 
and  if  M.  targets  or  more  are  present  after  the  (i  -  Oth^  attack,  then  M. 
missiles  are  f i red. 


CALCULATION  OF  THE  TRANSITION  KATR1-ES 

The  dependence  of  M. ,  p.,  b.  ,  and  T*  on  i  (the  attack  number)  will  be 
dropped  in  this  section  for  simplicity.  Let  the  elements  of  the  transition  matrix 
T  be: 

T. .  =  probability  that  exactly  i  targets  survive  given  that  j 
lJ  are  present  (;,j  =  Q,1,...n). 

For  purposes  of  use  in  Equation  (3),  it  will  be  understood  that  i  denotes  the 
column  and  j  the  row. 


Case  I:  Random  Targeting  with  Missiles  Spread  Among  All  Targets  Present 

This  is  the  most  difficult  transition  matrix  to  calculate.  We  have: 


T.  .  - 
'J 


0  i  f  j  -  i  <  0  or  j  -  i  >  M 

P. .  otherwise, 
i  J 


(9) 


where 


ptj  -  rM( 


i)Y — 1 

k  *  k  *  k  * 


ki +. . ,+k 


ki+1 


k. 


1  (1-s  1+1)...{l-s  J) 


;  do) 


s  =  1  -p  is  the  survival  probability,  and  the  sum  extends  over  all  ranging 

from  0  to  M  but  subject  to  the  constraint  that  their  sum  be  M(  i  k  = H) .  We 

r=l  r 

will  first  derive  this  result  then  simplify  it. 


Derivation  of  Transition  Probability 

We  are  looking  for  the  probability  that  i  targets  survive  if  j  are  pre¬ 
sent.  Label  the  j  targets  1  to  j.  Assume  that  targets  1  to  i  survive  after 

ki  ky 

being  targeted  by  k^  to  k.  missiles,  respectively,  (probability:  s  ...s  ) 
and  that  targets  i  +1  to  j  are  shot  down  after  being  targeted  by  k.+j  to  k. 


-k 


sum 


kj-s-i  k. 

mi  ss  i  1  es  ,  respect  i  ve  1  y  (probab  ility:  (l-s  )...(!  -s  ■').  The  k 1  s  must 
to  M  since  there  are  only  M  missiles  in  the  attack,,  There  are  M I  /  ( k  ^  ik^  .  „  .k.  i) 
different  ways  of  allocating  missiles  tc  target  1,  missiles  to  target  2, 

etc.,  and  kj  missiles  to  target  j.  Thus,  the  probability  of  a  given  com¬ 
bination  is  [Mi / ( k j Ikj I • • • k . ! ) ]  /  E  [ M i / ( k ^ i . k . I ] 9  which  is 
j  M! / (k . ! ! . . . kj ! )  according  to  the  multinomial  expansion.  Finally, 

there  are  ("|)  different  ways  of  choosing  the  i  survivors.  This  leads  to 
the  sum  in  equation  (10). 

Simplification  of  Transition  Probability 

The  sum  [equation  (10)]for  the  transition  probability  contains  an  exceedingly 
large  number  of  terms  and  must  be  simplified  for  actual  evaluation,  even  by 
computer.  The  basic  approach  in  simplifying  equation  (10)  is  to  get  the  summand 
to  reproduce  the  well-known  multinomial  expansion: 


(P1+p2+«--+Pj)  “ 


M 


•  k  1  k  . 

rK2'“°"  j 


kl  k2  kj 

:  Pi  P2  •••Pj 


where  I  means,  as  before,  the  sum  from  k  =0  to  k  =  M,  the  sum  of  the  k’s 

r  r 

always  being  M. 


The  quantity  we  want  to  reduce  is 


k  8  k  1  k  1 
,^2*  °  •  • 


;ki+-+k<o>h...o-skj> 


Consider  the  auxiliary  quantity 


F(b)  = 


k  1  U  1  k  1 

^1*^2* 900  j  * 


b .  k . 

. o ( 1  —  s  J  J),  ^ j) 


We  have,  of  course 


A  =  F  ( 1 ) 


03a) 


The  integral 


1 


1  -  sx  =  x  f  Ux-'dU 


(1M 


enables  us  to  put  F(b)  in  the  form 


F(b) 


/•••/  dui+l  •"  k  'k  ''  k  '  5  1 

s  S  A“" *  r  i .  k  2  - . . .  k  ■ . 


kn  +  . . ,+k  • 


X  b-  i  b.  k  •  ...  k.  u. 


i+1  **•  j  i+1  '**  i+1 


bi+lki+r1  ..bjkj"‘ 


...  u, 

J 


(13) 


where  sums  and  integrals  have  been  interchanged.  The  identity 

bk 


kubk  =  (Log  u)’1 


(16) 


puts  F(b)  in  the  form 


J  J  ,  .  . 

F(b)  dui+l  •**  duj  bi+l  •••  bj(ui+i  •••  Uj)  (Log  u)1_J 


—5^ —  _m: — 

3b.+-|.0.3bj  /Lji  k^!...kjl 


k1+...+ki  b.+]ki+1  b.k. 


ui+l 


(17) 
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The  sum  in.  equation  (17)  is  easily  recognized  as  the  multinomial  expansion  of 
b  :+i  b  •.  M 

(is  +  u.  ,  +  ...  +  u.  J)  .  The  derivatives  can  now  be  taken.  The  identity 

i  +  l  j 


M 


(a+u^)  *  M(a+uu)  uD  Log  u 


.b,^1  b 


(13) 


gives 


J-i  b.  ,  b. 

-  (Is  +  u.+1  +  ...  +  Uj  J) 


M 


3bT^i . . . ?b- 


(Log  u)J_1  M(M-1 ) . . . [M-(j-i)+l ] (is+u . , ,+1+. 
b.  ,  b.  1  ' 

uiU  ...  Uj  1  , 


b4i,  b. 


•+uj  J) 


(19) 


which  puts  F(b)  in  the  form 


1  J 


F(b)  =  M(M-l). . .  [M-(j-i)+l] /./ 


c*uj+1  •  •  -  bj 


b.,,-1  b  -r  s  b.+1  b. 

X  u-j+i  ***uj  ^ls  +  ui+l  +  ...  +  u  j  ) 


M-j+i 


(20) 


The  quantity  we  want  is  then 


.  ■  i 

A  =  F(1 )  =  M(M-1 ) . . .  [M-(  j-i  )+l  ]  (is  +  u.+1  +  ... 

M-j+i 


s  s 


+  Uj) 


dui+i  •••  duj  • 


(21) 
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The  integral  can  now  be  evaluated  b/  repeated  use  of 


1 


,¥ \n  =  (a+1  )pH'  -  (a+s 1 


dx(a-i-x) 


Ln+1 


n  +  1 


(22) 


The  result,  which  can  be  established  bv  induction,  is 


A  =  7  .  (-l)j-1"d  (Jdi)(d  +  (j-d)sj 
d=0 


M 


(23) 


The  transition  matrix  i^  chen'ore: 


T.  . 
U 


0  if  j  -  i  <  0  or  j  -  i  >  M 

J_M('j)  ^  (-1)J"1_d  (Jj1)ld+(j-d)s]  otherwise. 
d=0 


(24) 


This  form  for  T.^  is  far  more  tractable  than  equation  (9).  Indeed, 

equation  (10)  quickly  becomes  unwieldy  as  M  increases,  but  that  is  not  the 

case  for  equation  (23),  which  is  a  simple  sum  cf  not  more  than  n+  i  terms.  In 

other  words,  the  number  of  missiles  M  i s  no  longer  a  limit  on  the  "tractabi 1 i ty"  of 

the  problem.  The  number  of  targets  n  is  still  a  limitation,  however,  but  a  mild 

2 

one:  the  number  of  elements  in  the  transition  matrix  grows  as  n  , 


T  is  understood  to  be  1.  Note  that  along  the  diagonal  (i=j^0)T.. 

oo  n  1 J 

Note  also  that  Z  T.  .  =  1  as  expected.  Detailed  calculations  for  n=10 
i-n  1 J 

verified  that  equations  (9)  and  (24)  give  exactly  the  same  T... 


M 

=  s 
have 


Cas e  II:  Random  Targeting,  with  Missiles  Spread  Among  a  Specific  Number  of  Targets 
Let  us  now  assume  that  the  M  missiles  allocate  themselves  randomly  among 
only  the  closest  b  targets;.  Let  T  *  |  j  be  the  generic  transition  matrix  element 
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(probability  that  i  targets  survive  when  j  are  present).  If  b  >_  j  s  this 
restriction  makes  no  difference  and  there  is  no  change  in  the  transition  matrix 
T.  j  of  equation  (24) .  But  if  b  <  j ,  T’ .  ^  is  the  same  as  T.  _  ( j .  ju,  j  j  (j-b) 
Therefore,  the  full  transition  matrix  in  this  case  is: 


0  if  j  —  i  <  0  or  j  —  i  >  M 

T.  .  if  j  >b 
■J  - 

T 

i  +  b  -  j  ,  b  otherwise  , 


(25) 


with  T.  . 

i  J 


i  n  equation  (24)  . 


Case  IN:  Perfect  Targeting  —  One  Missile  Per  Target 

In  this  case  the  distribution  of  survivors  follows  a  simple  binomial 
distribution,  and  the  result  is: 


for 


( l”i  IP^O-P)"*3*1  for 

J  ’ 


for 


i  <_  j  <_  M 

M  <  j  <_  M+i  and  j^>  i 

M  +  1  <  j  or  j  <  i  . 


(26) 


For  example,  if  n  = 
arbitrary  o  is  (notation: 


10  and  M  =  5,  the  llxll 
a, be  means  ap^(l  —  p) c) : 


transition  matrix  for 


0 

0 

0 

0 

0 

0 

0 

0 

i-p 

0 

0 

0 

0 

0 

0 

0 

2,11 

(l~p)2 

0 

0 

* 

0 

0 

0 

3,21 

3,12 

(l-p>* 

0 

0 

0 

0 

4,31 

6,22 

4,13 

(1-K. 

0 

0 

0 

0 

5,41 

10,32 

10,23 

5,14 

(1-p)5 

0 

0 

0 

P5 

5,41 

10,32 

10,23 

5,14 

(1-p)5 

0 

0 

0 

p* 

5,41 

'10,32 

10,23 

5,14 

(1-P)* 

0 

0 

0 

P5 

5.41 

10,32 

10,23 

5,14 

(1-p) 

o 

0 

0 

P* 

5.41 

10,32 

10,23 

5,14 

0 

0 

0 

0 

p* 

5,41 

10.32 

10.23 

(27) 
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APL  programs  for  the  various  calculations  presented  here  are  available. 
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APPEND ! X  A 


APL  PROGRAMS  FOR  ATTRITION  CALCULATIONS 


APPENDIX  A 

APL  PROGRAMS  FOR  ATTRITION  CALCULATIONS 


This  appendix  contains  APL  programs,  written  by  Paul  E.  Klebe  and  the  author, 
which  calculate  the  attrition  based  on  the  method  described  in  the  text. 


The  first  listing  is  given  In  table  A-1.  Examples  for  each  of  the  three 
targeting  modes  are  given  in  tables  A-2  to  A-A.  In  the  case  of  random  allocation 
among  only  the  closest  b  targets,  only  the  case  b=M  has  been  coded.  Further, 
the  probability  of  kill  P  and  the  number  of  missiles  per  attack  M  are  the  same 
for  each  attack.  There  are  6  inputs  to  the  program: 

#  N  =  number  of  targets 

§  l<  =  number  of  attacks 

3  M  =  number  of  missiles  per  attack 

©  P  =  probability  a  missile  kills  its  target 

/  0  for  random  allocation  among  the  closest  M  targets 
'  1  for  random  allocation  among  all  targets  present 
2  for  perfect  targeting 


tTEST 


!1  if  N,  M,  P  and  TGT  have  not  changed  from  the  previous  run  (so  that 
transition  matrix  does  not  have  to  be  recalculated) 

0  otherwise. 


The  program  is  initiated  by  keying  in 
(TGT,  TEST)  RUNNOW  (N,K,M,P) 

with  the  appropriate  values  substituted  for  the  Variables.  If  it  is  desired 
to  make  it  interactive,  keying  in 


GETDATA 

will  cause  the  program  to  ask  for  the  required  inputs  by  name.  This  is  shown 
in  table  A-5. 


A- 1 


The  output  returned  by  the  program  is  as  follov<s: 

©  A  1 i st  of  i nputs 

©  The  transpose  of  the  transition  matrix  T  (not  T  itself  in  order  to  save  space) 
®  The  number  of  targets  surviving 
©  The  standard  deviation  of  targets  killed 
®  The  expected  number  of  missiles  fired 
®  The  initial  state  vector 
9  The  final  state  vector 

®  The  sum  of  final  state  vector  elements  (as  a  check:  it  should  be  1). 

This  program,  written  by  P.  Klebe,  is  an  original  improvement  over  an  earlier 
version  by  the  author  in  that,  in  order  to  save  space,  the  transition  matrix  T  is 
handled  as  a  vector  whose  elements  are  only  the  nonzero  elements  of  T.  This 
enables  the  program  to  run  for  large  values  of  N,  for  which  the  size  of  the  T 
matrix  is  prohibitive.  An  example  of  this  is  provided  in  table  A-6  ,  .he  the 
case  n=100  is  treated.  In  cases  when  N  is  i  arger  than  10,  only  the  n  ( zero 

elements  of  T  are  printed  in  the  output,  as  shown.  For  very  large  values  of 
N ,  the  number  limit  on  the  factorials  may  be  reached  and  the  program  will  not 
•-un.  Roundoff  errors  may  also  be  introduced  for  large  N,  due  to  the  very  small 
values  of  the  elements  of  the  final  state  vector.  (This  is  particularly  evident 
in  standard  deviation  calculations  )<, 

Running  time  is  of  the  order  of  seconds,  at  worst  minutes,  even  for  large 
N.  Note  that  under  perfect  targeting  and  random  targeting  among  closest  M„  no 
more  than  MK  targets  will  ever  get  targeted.  Thus,  if  ,  the  input  N  can 

te  chosen  equal  to  MK  and  N-MK  added  tc  the  number  of  survivors  given  by  the 
program.  This  reduces  the  run-time. 


Example:  N=50,  K=5 ,  M=6,  P=.l,  TGT=1.  The  program  correctly  calculates  ^7«3 
survivors,  but  the  final  state  vector  elements  sum  to  1.01,  which  causes  0 
in  (5)  to  be  -10.9.  This  is  why  a  statement  was  introduced  defining0*  to  be 
the  larger  of  the  calculated  value  and  0. 


A-2 


The  functions  used  in  this  program  are: 

®  The  overall  control  function  RUNNOW 
@  The  interactive  function  GETDATA 

#  TMR1  ,  which  calculates  the  transition  matrix  for  the  case  TGI  *=  0 

®  TMR2»  which  calculates  the  transition  matrix  for  the  case  TGT  53  ' 

o  TMR3»  which  calculates  the  transition  matrix  for  the  case  TGT  =  2 

@  SUM,  which  calculates  the  sum  in  equation  (24)  of  the  main  text 

©  NOMISP,  which  calculates  the  successive  state  vectors  and  number  of 
missiles  fired  under  perfect  targeting  assumptions 
©  NOMISR,  which  calculates  the  successive  state  vectors  and  number  of 
missiles  fired  under  random  targeting  assumptions 
9  BiGOATA,  which  handles  large  transitions  matrices  (N  larger  than  10) 

9  OUTPUT,  which  organizes  output  presentation. 

In  certain  cases,  however,  it  is  useful  to  be  able  to  handle  the  actus! 
transition  matrix  T,  as  opposed  to  a  vector  substitute  for  it.  The  following 
program,  written  by  the  author,  performs  the  same  functions  as  the  previous  one, 
but  handles  the  T  matrix  as  a  matrix.  Therefore,  it  cannot  handle  values  of 
N  that  are  larger  than  about  70  for  an  85,000-byte  workspace. 

There  are  4  inputs  to  the  program: 

9  N  =  number  of  targets 
©  K  =  number  of  attacks 
®  M  =  number  of  missiles  per  attack 
®  P  =  probability  a  missile  kills  its  target 
The  program  is  initiated  by  keying  in: 

PERTAR  for  perfect  targeting, 

RANTAR  for  random  targeting  among  all  targets  present,  and 
RANTARM  for  random  targeting  among  the  closest  M  targets. 

The  output  returned  by  the  program  is  the  same  as  before.  A  listing  is 
given  in  Table  A- 7  and  examples  for  the  3  targeting  modes  are  given  in  table 
A-8.  The  functions  used  are: 

©  The  3  overall  control  functions  PERTAR,  RANTAR,  RANTARM 
®  The  functions  calculating  the  state  vectors: 

NOMISP  for  perfect  targeting  and  NOMISR  for  re  .dom  targeting 
'9  The  3  functions  TMi,  TM2 ,  and  TM3,  calculating  the  transition  matrices 
for  the  3  targeting  modes:  random  among  closest  M,  random  among  all,  and 
perfect,  respectively. 

A- 3 


TABLE  A- l 


LISTING  QF  FIRST  Af’L  PROGRAM 


) LOAD  KLEBE 

SAVED  B0/10/2P  17.00.13 
)FNS 

BIODATA  GETDATA  NOMISP  NOMISR  OUTPUT  RUNNOW  SUM  TMR1  TMR2  TMR3 


VBIGDATACD3V 
?  BIGDATA  »  KC  *  CT . VT 
C  1  3  SPACE 

C23  1  VALUES  OF  N.M.SIZE  OF  T' 

C3  3  N.M.ipT 

F4  3  SPACE 

FS3  >(TEST=1 )pO 
FAT  CAi-KCl-0 

F73  ' THE  TRANSPOSE  OF  THE  TRANSITION  MATRIX  T  ISt  1 

FBI  FIERI  IN  8  VT!  H  Ml. N-CA 

E93  '  IAJF7.4' ?<CA»TFKC+WT3> 

C  103  KC*  KC+VT 
Fill  ->  (  NiCA<  CAT  1 )  p  BERLIN 
v 

^GETDATAFOlv 
v  GETDATA 
FI  3  SPACED'  ' 

C23  0I0<-0 

F.33  'LIST  THE  NUMBER  OF  TARGETS.  NUMBER  OF  SUCCESSIVE  ATTACKS.  ' 

F 4 3  '  THE  NUMBER  OF  MISSILES  IN  EACH  ATTACK.  THE  MISSILE  PROBABILITY  OF  KILL' 

T53  F'EST«-a 
CA3  SPACE 

F73  'TYPE  1  FOR  PERFECT  TARGETING.  2  FOR  RANDOM  TARGETING’ 

TOD  TGTf-3-0 

F93  ->  <  TGT ”2 >  p BOSTON 

C 1 03  SPACF 

FI  13  'TYPE  1  IF  MISSILES  SPREAD  AMONG  ALL  AVAILADLE  TARGETS' 

F.  123  '  2  IP  MISSILES  SPREAD  AMONG  A  SUBGROUP  OF  M  TARGETS’ 

F  1 33  TGT «-2-  0 

F 1 4 3  BOSTON: SPACE 
C 1  S3  SPACE 

F163  'TYPE  1  IE  THE  T  MATRIX  IS  SAME  AS  IN  PREVIOUS  RUN  —  OTHERWISE  TYPE  0’ 

C  1 73  TEST<-0 

F103  (TGT. TEST >RUNNOW  PEST 
v 


vNOMlSPCm? 
v  NOMISP. AT? VT.JJ 

m  Em-o 

F23  ATH 
C33  ST<-SO 
F43  MN<-MLN 
F53  GG*RT>1- i/STF >MN3 
CA3  FM< EMT(MNxRT)+f/STC\MN3*\MN 
F7  3  JJ<KT>0 
FBI  ROME J VT< 1TMLN-JJ 
[93  STEJJ3^T/STCJJ+\VTDxTCKC+\VT3 
C103  KC^KCTVT 
C113  -*  (Ns  JJ«- JJ+1  )*>ROME 

C  123  -»<KiAT<  AT+1  >*GG 
v 


A-  k 


A-l  (con't) 


*  NOMISR?AT?VTSJJ 
Cl  3  EM*© 

122  ATM 
m  «rr*-sn 

C43  6G?EM<-EMTMxl~STC©3 

C53  JJ^KCfrO 

r.62  ROMES  VTM  EMIN- JJ 

C73  STEJJ3M/STEJJ4-WT3XTCKC+IVT3 

l83  KOKC4VT 

C93  CN'i  JJ<- JJ+1 )  pROME 

C103  -»<KS:AT*AT+1  >pGG 

^OUTPUTCOl^ 

»  OUTPUT 

Ell  NAME* 3  32 p  '  RANDOM  TARGETING  AMONG  M  CLQSESTRANDOM  TARGETING  AMONG  ALL 

C23  SPACE 

C33  NAMECTGTM  1  PERFECT  TARGETING 

£4  3  SPACE  L» 

C53  ' THE  INITIAL  NUMBER  OF  TARGFTS  IS  'rN? 

C63  ® THE  NUMBER  OF  SUCCESSIVE  ATTACKS  IS  '  ?K? 

C73  'THE  NUMBER  OF  MISSILES  IN  EACH  ATTACK  IS  ' f M? 

C81  'EACH  MISSILE  KILLS  ITS  TARGET  WITH  PROBABILITY  *  IP? 
t?2  SPACE 

EI03  -*<  <NrM)>10)pPARIS 
CIO  KC<CA<-0 

023  'THE  TRANSPOSE  OF  THE  TRANSITION  MATRIX  T  ISI  e 
033  SPACE 

043  MADRID?  VTM  +  MLN-CA 
053  KUNTl-CAfVT 

063  9  5? (CApO ) ? TCKCT * VT3 »  <OfKT) pO 

071  KC*KCFVT 
003  -KNsCA^CA+DpMABRID 
093  PARIS 5 SPACE 

E203  ‘THE  NUMBER  OF  TARGETS  SURVIVING  IS  '  ?E? 

120  'THE  STANDARD  DEVIATION  OF  TARGETS  KILLED  IS  ' PSD! 

C223  'THE  EXPECTED  NUMBER  OF  MISSILES  FIRED  IS  1 fEMf 

C233  SPACE 

C213  'THE  INI I T I AL  STATE  IS  ' 

C253  1©  7?S0 

C261  SPACE 

L273  'THE  FINAL  STATE  IS  ' 

C231  1©  7?ST 

C293  SPACE 

C303  ' A1 JF6.2' t ( 'SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  —  IT  IS?  ' f+/ST> 

C313  -K<NrM)«ll)pO 
C323  BIGDATA 
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A- 1  (con't) 


vRUNNOUCOlv 

V  BAKER  RUNN0W  ABLE*ST2 
0  3  TGT*-BAKERC01 
E23  TEST  <-B  AKERE  1 3 

E33  MAT1+-3  4p 1 TMR1TMR2TMR3 ’ 

E43  N4-ABLEC03 
E53  K<- ABLE  Cl  3 
E&3  M<~ABLEC23 
E73  P<  ABLEC33 

C81  MAT24-3  6p  1  NQiIISRNQNISRNOMISP  * 

C91  SO<- (NpO  )  » 1 

003  -MTEST=1)pL0ND0N 
013  ®MAT1.rT3Tn 
023  LONDON  i  SPACE 
033  SMAT2ETGT  >  3 
043  ]>iN+l 
053  E<-+/IxST 

063  ST2«-<-f/<T.*2)xST)-(-f/XxST)*2 

0  73  ST2f-ST2rO 
083  SB<-ST2*0«5 

093  OUTPUT 

v 

vSUMCnjv 

v  SIG<I  SUM  JJVARrCl  iBlf  B 
03  DHUCltJ-I 
C23  Dl«-D-fJ 
C33  VARf-~l+2x2  I  Cl  +  t 

C43  SIG4-  ( I  ♦  J )  X-/VARX  ( D« Cl  >  x (Dl+Sx 1~D1 ) ftfl 
v 

^TMRl  CD3'? 
v  TMR1 ?S?I»J{V 

03  T«-(L0.3FN+H(Mx0rN"M>  +  (  ( 1+S)  xSfrMLNH-2)  pO 

E23 

r.33  SK4-<S<-1-P)*M 

C43  L.OOPJJM  +  1 

C53  TEV3*SK 

C63  V<-V+l 

E73  -K  I=N) pSTEP 

E83  LOPJ-KJ>I+M)pSTEP 

E93  TCV3<-(I-0r J-M) SUM  JLM 

003  V<-V-H 

013  >  <  Ni  J<~  J+l  )  pLOP 

023  STEPJ-KNOt-I  +  DpLOQP 

033  TC03O 
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A- 1  (con 1 1) 


^THR2CD1^ 

V  TMR2S SI  1 » J» V 

E13  T«-(L0.3+N+i  +  <MX0rN-M>  +  (  <1+S>xS<-MLN)*2)fO 

t2z 

E33  SK*-  (S<-l-P )  *M 

E43  LOOP*  J<~  1  +  1 

C53  TCV3<-SK 

C63  V4-V+1 

£71  -><I=N)pST£P 

C83  LOPS-><J>I+M)pSTEP 

C93  TCV3+-I  SUM  J 

C103  V«-V+l 

CU]  -»<N*J<-J+l>pLaP 

C 123  STEP?'KNsI<-I  +  1)pLOQP 

C133  TCQ3M 

V 

?TMR3C£]3v 
^  TMR3p3» I » J» V 

Cl  3  T<-  <  L©«  3+N+1+  (MxOrN-M  >+<(t+S)  XS<-MLN  )->2 )  pO 

C23  V<-I<-0 

C33  'SK<'<Sfl--P)*li 

E43  LOOP*  JM 

C53  LOP*-* <  J> I+M) pSTEP 

C63  TCV3<  ( P*J~I ) x  < ( J-I > i JLM ) xS*I-0r  J-M 
C73  V«-V+l 
C83  -><NsJ*-J+1)pLOP 
C93  STEP1-KNj:I*I  +  1>pL0QP 
C103  TC03M 
v 
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TABLE  A- 2 


EXAMPLE  OF  FIRST  PROGRAM  FOR  RANDOM 
TARGETING  AMONG  M  CLOSEST  TARGETS 


(0,0)  runnow  uo>4>5'.3)  key  In 


RANDOM  TARGETING  AMONG  M  CLOSE sV 

THE  INITIAL  NUMBER  OF  TARGETS  IS  1© 

THE  NUMBER  Of  SUCCESSIVE  ATTACKS  IS  4 

THE  NUMBER  OF  MISSII  FS  IN  EACH  ATTACK  IS  5 

EACH  MISSILE  KILLS  ITS  TARGET  WITH  PROBABILITY  0.5 


THE  TRANSPOSE  OF  THE  TRANSITION  MATRIX  T  IS] 


1.00000 

0 .96075 

0.55664 

0. 1501 R 

0.02197 

0.00120 

0.00000 

0.00000 

6 . OOOOO 

0.00000 

0.00000 

0.00000 

0.03125 

0.41211 

0 . 50926 

0.22339 

0.04200 

0.00120 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.30131 

0.46692 

0.2E500 

0.04200 

0.00120 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03123 

0.25647 

0.43000 

0.25500 

0.04200 

0.00120 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.23255 

0.43000 

0.25500 

0.04200 

0. 00120 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.03125 

0.23255 

0.438AO 

0.23500 

0.04200 

0.00120 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.23255 

O.43R00 

0.25500 

0.04200 

0.00000 

0.00000 

0.00000 

0.00000 

0 . OOOOO 

0.00000 

0.00000 

0.03125 

0.23255 

0.43800 

0.25500 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.23255 

0.43800 

0,00000 

0.00000 

0.00000 

0.00000 

0 . OOOOO 

0.00000 

0.00000 

O.i-'OOOO 

C. OOOOO 

0.03125 

0.23255 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . OOOOO 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

THr  NUMBER  OF  TARGETS  SURVIVING  IS  2.04 

THE  STANDARD  DEVIATION  OF  TARGETS  KILLED  IS  1.49 

THE  FXPECTED  NUMBER  OF  MISSILFS  FIRED  IS  20 

THE  INJITIAL  STATE  IS 

0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0090009  0.0000090  0.0000090  0.0000000  1.0000000 
THr  FINAL  STATE  IS 

0.1654055  0,2340347  0. 2455501  0,3863430  0.1065119  0.0453583  0.0136313  0.0027969  0.0003703  0.0000284  0.0000010 
SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  —  IT  IS]  1.00 
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TABLE  A- 3 


EXAMPLE  OF  FIRST  PROGRAM  FOR  RANDOM 
TARGETING  AMONG  ALL  TARGETS  PRESENT 


<i,o>  fhjwwou  (i©.4»5?.s>  key  in 


RANDOM  TARGETING  AMONG  ALL 

THE  INITIAL  NUMBER  OF  TARGETS  IS  10 

THE  NIIMRFR  OF  SUCCESSIVE  ATTACKS  IS  4 

IMF  NUMBER  OF  MISSILES  IN  EACH  ATTACK  IS  S 

EACH  MISSILE  KILLS  ITS  TARGET  WITH  PROBABILITY  0.5 


THE  TRANSPOSE  OF  THE  TRANSITION  MATRIX  T  IS! 


1 .00000 

0.96075 

0.55664 

0. 15010 

0.07197 

0.00120 

0.00000 

0.00000 

0.00000 

0,00000 

0 

ooooo 

0.00000 

0.03125 

0.41211 

0.50926 

0.22339 

0.04200 

0.00289 

0. ooooo 

0.00000 

0,00000 

0 

ooooo 

0.00000 

0.00000 

0.03125 

0.30131 

0.46692 

0.25500 

0.05707 

0.00469 

0.00000 

0.00000 

0 

ooooo 

0 . 00000 

0.00000 

0.00000 

0.03125 

0.25647 

0 . 43000 

0.27247 

0.07028 

0.00641 

0.00000 

0 

0000© 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0 . 23255 

0.41775 

0.20309 

0.03011 

0.00800 

0 

ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0. OOOOO 

0.03125 

0.21776 

0.40296 

0.20999 

0.0000? 

0 

00945 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.20774 

0.39173 

0.29473 

0 

09450 

0.00000 

0.00000 

0.00000 

0,00000 

0.00000 

0.00000 

0. OOOOO 

0.03125 

0.20051 

0.30294 

0 

29013 

0.00000 

0.00000 

0 , ooooo 

0.00000 

0. ooooo 

0.00000 

0.00000 

0.00000 

0.03125 

0. I 9505 

0 

37589 

0.00000 

0.00000 

0.00000 

0.00000 

0,00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0 

19070 

0.00000 

0.00000 

0.00000 

0.  0000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 

03125 

THE  NUMBER  OF  TARGETS  SURVIVING  IS  1.72 

THE  STANDARD  DEVIATION  OF  TARGETS  KILLED  XS  1.44 

THE  EXPECTED  NUMBER  OF  MISSILES  FIRED  IS  19.9 


THE  INTITJAL  STATE  IS 

0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0600000  0.0000000  1.0000000 
THE  FINAL  STATE  IS 

0.73S0293  Q.257S490  0.2312340  0.1553414  0.0790975  0.0305000  0.0097932  0.0018388  0.0002639  0.0000233  0.0000010 
SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  ~  IT  IS!  1.00 


TABLE  A” 4 


EXAMPLE  OF  FIRST  PROGRAM  FOR  PERFECT  TARGETING 


(2,0)  RUNNDU  UQr4»5».5) 


PERFECT  TARGETING  RA 

THE  INITIAL  NUMBFR  OF  TARGETS  IS  10 

THr  NUMBER  OF  SUCCESSIVE  ATTACKS  IS  4 

THE  NUMBER  OF  MISSILES  IN  EACH  ATTACK  IS  5 

EACH  MISSILE  KILLS  ITS  TARGET  WITH  PROBABILITY  0.5 

THE  TRANSPOSE  CF  THE  TRANSITION  MATRIX  T  IS! 


1 .00000 

0.50000 

0.25000 

0. 12500 

0.06250 

0.03125 

0.00800 

0,00800 

0,00000 

0.00000 

0.00000 

0.00000 

0.50000 

0.50000 

0.37500 

0.25000 

0.15625 

0.03125 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 25000 

0 . 37500 

0.37500 

0.31250 

0.15625 

0.03125 

0.00000 

0.00000 

0.00000 

0  00000 

0.00000 

0.00000 

0.12500 

0.2S000 

0 . 31 250 

0.31250 

0.15625 

0.03125 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

o.ooooo 

0.06250 

0.15625 

0.31250 

0.31250 

0.15625 

0.03125 

0,00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.15625 

0.31250 

0.31250 

0.15625 

0.03125 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.15625 

0.31250 

0.31250 

0.15625 

0.00000 

o.coooo 

0.00000 

O.OOOuO 

0.00000 

0.00000 

0.00000 

0.03125 

0.15625 

0.31250 

0.31250 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.15625 

0.31250 

0 . 00000 

0.00000 

0. 00000 

0.00000 

o.ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.03125 

0.15625 

0.00000 

0.00000 

0.00000 

6.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.03125 

THE  NUMBER  OF  TARGETS  SURVIVING  IS  1.44 

THE  PI  ANBARB  IlEVIATinH  OF  TARGETS  KILLED  IS  1.23 

THF  EXPECTED  NUMBER  OF  MISSILES  FIRED  IS  17.1 


THE  INIXTIAL  STATE  IS 

0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0,0000000  1.0000000 

THE  FINAL  STATE  IS 

0.2377130  0.3460693  0.2413464  0.1120377  0.0421301  0.0147858  0.0046206  0.0010G72  0.0001812  0.0000191  0.0000010 
SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  —  IT  IS!  1.00 
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TABLE  A- 5 


EXAMPLE  OF  INTERACTIVE  VARIANT  OF  FIRST  APL  PROGRAM3 


©►  GETOATA 

LIST  THE  NUMBER  OF  TARGETS.  NUMBER  OF  SUCCESSIVE  ATTACKS. 

THE  NUMBER  OF  MISSILES  IN  EACH  ATTACK.  THE  MISSILE  PROBABILITY  OF  KILL 

0.* 

^  10.4. E. .5 

TYPE  1  FOR  PERFECT  TARGETING.  2  FOR  RANDOM  TARGETING 

n: 

2 

TYPE  1  IF  MISSILES  SPREAD  AMONG  ALL  AVAILABLE  TARGETS 

2  IF  MISSILES  SPREAD  AMONG  A  SUBGROUP  OF  M  TARGETS 

□  : 

ss-  2 


TYPE  1  IF  THE  T  MATRIX  IS  SAME  AS  IN  PREVIOUS  RUN  —  OTHERWISE  TYPE  0 

n: 

t*-  o 


RANDOM  TARGETING  AMONG  M  CLOSEST s etc .  (Output  Is  as  in  Table  A-2). 


^  means  "key  in". 


A-l  1 


TABLE  A-6 


EXAMPLE  OF  FIRST  PROGRAM  FOR  LARGE  NUMBER 
OF  TARGETS  (N  >  iO)s 


(0-»J>  nilNMOU  (  100»  >0i5i  .  I )  Key  I  n 


|M  I  Af'f'J  1  I  fJF»  M  Cl  Gr»f  ST 

n-  u-'Mr.'it  ni  .h  i  u  nr  TAt,r:rrr  rr  100 

t*  r  tr  cm  rnrrr?.r.iyF  -m  tacks  is  ?0 

'mt  rn«y:rj  in  Mir.r  u  1 1;  tn  tami  attack  is  5 

r.  r-kj  rur  K?jt  r.  nr>  taugf t  with  rftfWAftn  nr  0.1 

t»  t  ( r.  nr  rnr.r.rfn  '•"hvivinr  is  ?^.-y 

t  rr  r>  u  y  at  ?  ah  or  Terris  mi  i  r.n  is  2.03 

ir*  1  yi  1  nt.ri  m'-mitk  01  Mr.j.nrr.  rirn»  is  100 

If!  T  r*  I  I  T  T  ft!  STMT  I  r. 

o.sA^.r.ooo  n  r  ^or.oo r..  *00*000  c  .0000000  0,00000 oo  o. 0000*00  0,  .00000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000 

.'.■'W  n.'-r.-vr.-.o  0.000f0?0  0,0000000  0 . 0"*o*00  0 *0*00000  a,<-.o0000O  0  OOOOCOO  O.OOOCOOO  0. (>000000  0.0000000  o.oo 

OOAOO  A  -.v.r  A  ,  Or,  00  0  10  0, 000000 C  O.fOOOOOO  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000 

■  > .  OOOC-A  v*  ii/,nr,.ir.GA  O. 0000000  O  ,  <*  <  ’0*0*0  O.  0»:/'0Cl/'/'  0 . 00A<  (u  C  0.0000000  0.00*00*0  0.00*0000  0.0000000  O.OOoOOOO  0.00 
O'.'1'OO  I.."  ■  ■' .  .•0'-'i;-000  (  .  O  ,.\y  •’■<»*  o,  *n<V)o.ip  O .  0*-V#  >*•'  o. f 000*00  0.  -.>000000  0.0*00000  0.0000000  0.0000000  0  ,  *00*000 

o.rorr-i  on  0.0  i  A  .’i  r*.  .’I  (  ,  (*OO0u0O  0  .  **0C  **0  0. 00*000 j  O.OOOOO'lO  0. (‘*00000  0.000*00*  o.ooooooo  0.00000*0  0.0000000  0,00 

n.n'ir.o  .  O  .  *i\‘  o.Oi-O  O  ,  0000  'TV'.  *.'0**0*.'  0. 800*000  O  .  ****  (v«0  O.0O00000  *.0000  AOO  0.0000000  0.0000000  0.0000000 

•J  .  0  A  o  .  •:  ’'.'.O'V  •>  *.000**0*  0.0000 ’-00  0.0«vooon  0.0000000  *, **00000  0.0000000  0.0000000  0.0000000  0.0000000  0.00 

00**0  0,t.-*0000*  o  .  ."AOOOOO  *..H‘00>*0  0  ,  OOOOOt- '  0  .  **■'•'  u>00  0.0r.“00*0  0.0000000  l.OOOOOUO 

ns  f  T»;oi  riAir  ir. 

on  0 .  0000***  n. *0**000  O.OOAO**')  A  ,0000000  0.l**v'000  0.0000000  0.0000*00  0.0000000  0.0000000  0,0**0000  0 .  OOOOCOO 
.  Of.rc.Ot'O  f',-\oT00.)  0.000'0,'T  *  o»W00*  *,0*0<?  ;**•*  O.OOOOOOO  O.OOOOOOO  0.0000000  *.**00000  0.0000000  0.0000000  0*00 
A(.  I.  nA.iA  v*  .  o.  >AAA(>00  0.  WW.*-"*  .  <‘0****0  0 . 0(m-0- >  A .  CO'W.' *  0.0000000  A* 0*0*0**  0 . 0*00000  0.0000000  0.0000000 

•■..iVii.  ,  r(.*'v.,>A',  A. ooo*  «  M'  ,.fin,.Oo  *.*'•  «■>.»' n.x-OOO  “*.*0*0*00  0.0000000  0.0000000  *.**0**0©  0.0000000  A. 00 

A.m.aa  ,  s.VO-  •  '■  .  OuO  '  OOA  o.r  •***•<  O.C.'OOO.'-'  o .  oo^'.vko  o.ooooooo  O.OOOOOOO  0.0000000  O.OOOOOAO  o.ooooooo 

o,  •  .’»/.  "»  .;.(V0'"  •'  o.rt.vn  O.o.  »r.v.i  ••.00^0000  O  or.OoOOo  0.00i)*'‘>00  O.OOOOOOO  O.OOOOi'A  0.0000000  0.00 

c  V  AO  0  .  .H'O0< 1  O.OT-  .'OCO  .r-.  0  .  o.-.-OO-’l  0.0>OA.-.oi  o.oooonj*,  O.OOOOOS?.  o.  00001  f:0  0.00(»0!i  T*  O.OOOi6?t  O.OOO'.AAi 

A  O.i  J  J  •••  •:  O.O'.IVI  o.oo/.'  *  I  o.o*  *»  •/,;«•*  O .  or*  O’Oi*  *  O.iMJf'.'i  O.  M  O-lSUO'i  O.M0I?r»?  0.13 

71  1  "i!  0.’0A*‘T*  O.O/./V/M  O.OtVO/49  O .  i*  0/.*.* .  :  ••  0 . 001  '•*!  .  0.00070  ISS  0.0'003/./» 

r”M  nr  i  ;m/>i  r.TATi  i:  in  nit;  nunm  i.  j:»  i  ri  ir:  t.^o 


f;ote:  only  the  nonzero  elements  of  the  transition  matrix  are  printed.  To 
explain  their  location,  it  is  best  to  give  an  example.  Below  is  how  the  T 
matrix  in  table  A-2  would  be  printed  in  this  format  (obtained  by  keying  in 
BIGDATA  after  program  execution). 

BIODATA 


VALUES  OF  N» SIZE  OF  T 
*0  5  51 


TRANSPOSE 

OF  THE 

TRANSITION  MATRIX  T 

is: 

0 

1.0000 

0.9688 

0  *  5566 

0,1582 

0,0220 

0,001.2 

i 

0,0313 

0.4121 

0.5093 

0,2234 

0.0420 

0,0012 

2 

0.0313 

0.3013 

0,4669 

0,2550 

0,0420 

0,0012 

3 

••V. 

r-t 

a 

o 

0 , 2565 

0,4300 

0 , 2550 

0.0420 

0,0012 

4 

0,0313 

0,2326 

0.4380 

0 . 2550 

0,0420 

0,0012 

5 

0,0313 

0,2326 

0,4380 

0 , 2550 

0,0420 

0,0012 

6 

0,0313 

0,2326 

0 , 4380 

0.2550 

0,0420 

7 

0,031.3 

0,2326 

0,4380 

0,2550 

8 

0,0313 

0,2326 

0,4380 

9 

0.0313 

0,2326 

10  0.0313 


AI2 


Values  of  N.M,  'Size  of  T 
100  5  591 

1  *  •  !  .  •  II  ■  l  !  t:-\  *  I  '  •  *•  : 


1 

4  • .  • 

•  • .  1 

o.--'  -.) 

f:  ,  •  ' 

0.*:  •■'! 

<:  f  .. 

<•  .  •  '  * 

•  .  \  * 

f,.'»  ••■>  1 

.1  .  -.A 

4 

0  .  •«  *. 

O .  •  *7 

1 

•*-.  ■  ?0 

■  .  v. •*■■ .  <.«,**'•■ 

0.*  •< 

•  .  '  1 

6  50 

o. i*.  :  *  ••• 

'  • 

•  ■  .  1 

O..V.1 1 

O.IM-O 

*'  .  <  ■ .*  ; 

c.  .  ■  t  ' 

1 .  '  1 

.  An.  0 

•  6  .  -  V 

?>.'/•  *■ 

.  • .-  :  • 

O.I, (  -1 

U.O'-*'f> 

.  *.W  .  6  .  M 

0.*  ' 

•  ;17 

o.t  (  f«1 

0,0600 

JO 

o,'.,vo'>  ».:v.  •*.) 

0 . 063  * 

0.0047 

0.000! 

0 . 0000 

1  1 

o.1  'o 

0.0673 

0.  C/'V 

A  # 

0.0000 

1  7 

o.*  V-  •;  o.  »4.',o 

0.0' 77 

0.00  17 

0.000*. 

<V  .  0900 

1  ' 

o.*.*  0.  '4  ’0 

0.  7 7 

0,0047 

0.000! 

0.0000 

!  1 

o,*  '• .  '  1  7  A 

0  ». 

A.  0001 

0 . 00 AO 

1 

o. n,.jno 

0.063  r, 

0,6047 

0  .  AOf'l 

0.0000 

1/ 

6.5765  0.14  70 

0.0673 

0,0047 

0.0001 

0.0000 

1  / 

O.‘vc.0  0.14  70 

0.0673 

0,4417 

0.0001 

0.0000 

1  !' 

A.  5  0.74  7  0 

0.4'.?.! 

A .  '‘•A47 

O.OoOl 

0.0000 

1  O 

0.7470 

0.06*»3 

0,0047 

A . AAO J 

0.0000 

70 

O.'VV.  0.74  7  A 

0 . 067.7 

0 , 0A4O 

0.0001 

0.0000 

:  i 

0 . VO*'  0.74  70 

0.0637 

0.004? 

0.0001 

0.0000 

:*.» 

0.5905  0.7170 

0,067.7 

0.0047 

0.0001 

0.0000 

■.  *. 

n.f.w,  0.74  70 

0.0677. 

O.O04  7 

0.0001 

0.000" 

O.f.r^r.  A.  'VO 

0.0673 

0.004? 

0.0001 

o.oroo 

?5 

0  .  5905  0  •  ,7470 

0.0673 

0 • 004? 

A. 0001 

0.0000 

:v, 

0.5905  0.7470 

0.0671 

0.004? 

o.ooot 

0.0000 

:*7 

O.!.V0t>  0.34  70 

0.0673 

0.0042 

0.0001 

0.0000 

O.rn-jr.  0.34*0 

0.0672 

0.O04'1 

0.0001 

0.0000 

:  ? 

A.57C5  0.  »  1.70 

a  , .  • .( 

o.A.vtr 

0.0001 

0.0000 

■0 

0 .5.905  0.7470 

0.0623 

0.0047 

A. A001 

0.0000 

» • 

Q.591-r'  0.3170 

0 . 0673 

0.AA4? 

0.000* 

0.0000 

7.  "J 

0.7.707.  0,3 4 7r) 

0.0673 

0.004? 

0.0001 

0.0000 

i  / 

0.5705  0.34.70 

0.0677 

0.004? 

0.0001 

0,0000 

7-1 

o.5?cr;  o.n^o 

0.0673 

0.004? 

c .0001 

0.0000 

35 

o.r.9on  0.3430 

a. 0633 

0 . 004? 

0.0001 

0.0000 

30 

0.590B  0.7430 

0.0672 

0.004? 

0.0001 

0.0000 

37 

0.5*05  0.24  30 

0.0673 

0.0042 

0.0001 

0,0000 

35 

0.5‘05  0.3470 

0 . 0623 

0.0042 

0.0001 

0.0000 

39 

G.V>05  0.34  70 

0.0633 

0.0042 

0.0001 

0.0000 

40 

0.!I90S  0.74  70 

0.0632 

0.004? 

0.0001 

0.0000 

4  1 

0.3905  0,34.30 

0.0637 

0.00-12 

0.0001 

0.0000 

/  ? 

0.V?05  0.34  30 

0.0673 

0.0047 

0.0001 

0.0000 

4  < 

o.'ivrj  o.7vo 

0.0*73 

0.0042 

0.0001 

0.0000 

44 

0. •;°o*)  0.34.70 

0.0633 

0.004? 

0,0001 

0.0000 

4*5 

0 . 5°0!i  0.J470 

0.0673 

0.004? 

A . 000 1 

0.0000 

46 

0.590?  0.347C 

0.0671 

0 .  OO.i  •» 

0.0*01 

0.0000 

4  7 

0.5*>0ri  0.3470 

C.  0.633 

0.004? 

0.0001 

0.0000 

.IP 

0.5*0*“  0.7470 

0.0633 

O.0/' 

A  000 1 

0.0060 

4? 

0.r.90*l  O.’.M’O 

0 .0^  72 

0.0042 

0.0001 

0.0000 

r.o 

0.7:90S  0.34.10 

0.06"*  3 

0.0042 

0.0001 

0.0000 

r  1 

A.V7 or.  0.747.0 

0,0677 

0.C-A4? 

0.0001 

0.0000 

r.? 

0,5?**  5  0.2470 

0.0673 

0.0042 

0.0001 

0.0000 

0.590:.  C . 7 4 7. 0 

0 . 0637 

O.C042 

0.0001 

0.0000 

f,.? 

0.5705  0.2.420 

0.0473 

0.0042 

0.0001 

0.0000 

r.j 

0.5905  \J . 3430 

0.0633 

0.C04? 

0,0001 

0.0000 

:»i 

0.3705  0.34.70 

0,0677. 

0.0042 

0.0001 

0.0000 

f,  7 

o.s*- or.  0.7.430 

0.0633 

0.0042 

0.0001 

0.0000 

!'fS 

0  5905  0.3470 

0,0673 

0.004° 

0.0001 

0 . 0000 

B-y 

0.3905  0.3470 

0.0633 

0.0042 

0.0001 

0.0000 

/0 

0.5905  0.347.0 

0.0633 

0.0042 

0 . 000 1 

0.0000 

At 

0.5*05  0.. 7-1/0 

0.0633 

0.0-142 

0,0001 

0 . 0000 

/.? 

0.59 Dr.  0 . 3 4 30 

0,0673 

0.0047 

0.0001 

0.0000 

A7 

0.9405  0.7*70 

0.0673 

0,0042 

0.0001 

0  .0000 

64 

0.7.905  0.3430 

0.0673 

0.0042 

0.0001 

0.0000 

AU 

0.5905  0.7470 

00/.**  3 

0.004? 

0.0001 

0.0000 

A  A 

0.5905  0.74.79 

0.0623 

0.0042 

0.0001 

0.0000 

67 

0.5°<:5  0.34*»0 

0.0633 

0,004? 

A. 0001 

0.0000 

60 

0.5905  0.7.130 

O,067*» 

0.0017 

0.0001 

0.0000 

69 

0.5905  0. 3**30 

0, 0673 

0.004? 

0.0001 

0.0000 

70 

0.5 V05  0.3430 

0.0673 

0.00*? 

o.ooci 

0.0000 

71 

0.5905  0.3470 

0.0677 

0.0C4? 

0.0001 

0.0000 

77 

o.r°or.  0.7.'  *o 

0 . 0673 

0.0042 

0.0001 

0 . 0000 

7.3 

0.5905  0.2470 

0.0677. 

0  .OA-17 

0  0001 

0 . 0000 

74 

0.5905  0.3430 

0.0677. 

0,0042 

0.0001 

0 . 0000 

7C 

0.5905  0.  1430 

0.0673 

0.004? 

0.0001 

0.0000 

0,5'': 5  0.34  30 

0.0673 

0.004? 

0.0001 

0,0000 

77 

o.5?05  0.3430 

0,0673 

0 .  r-M? 

0.0001 

0.0000 

70 

A.r.yo,',  o.343o 

r> .  0/  73 

o.oo-i;? 

o.ooot 

0.0000 

79 

O.r.^O*.  0.347.  •• 

0.0/77. 

0.06-1** 

0.0601 

0.0000 

110 

0.56*05  0.3430 

0.0673 

O.O.'I? 

r ,0001 

A  .  ,1000 

111 

0.5905  0.34:«<> 

0.0/.73 

fi.cr 

A.  (>AA| 

0.0000 

n? 

0.5° 05  0.  *.4  30 

0,0633 

0.00 1? 

O.OOt  J 

0.0000 

B3 

0.5905  0.3HO 

O.a- 73 

0.0  74  7 

A , AA^ t 

0 . 0006 

f'4 

o.r.vy  0..7470 

0,0.673 

0.0017 

0 .000 , 

0.0000 

Oli 

0.5 70 5  0,31.5' 

0 , 0.-.72 

0 .  fi..  1  ? 

o.ooot 

0.0000 

i V. 

0.1.905  0.3.4 .!t' 

0. 0.673 

0 .  O’'-’? 

0,0001 

0.0000 

l!  .* 

o.r.sv'.  0 .  51  *0 

A.  or  1  7 

0.0001 

r  .co.-o 

I 

■'.590*  0.34  7.0 

O.Oa77 

A, 0047 

,»  A,>f|  , 

-...'OOO 

V" 

O.W.  0,31.5} 

O.a.^I 

o.o.'l- 

0.  ‘  )0I 

A. AO AO 

•■*  O 

o.:.Vi*5  0., 

0,0  ? 

0 .  1? 

C.  t'OOl 

.1 . 0  A OO 

0.5  V5  0.  '-1  '  \ 

0,0. 73 

O.OA-17 

0 .  orr.  1 

A  ,  i'Pi'O 

w.O  .7  f. 

4.0.'  1/ 

0  .  <  1 

0.-'*'*V 

'•  * 

'".59 ,7*.  0.“S*  ■ 

A , ; 

<• ,  *■  *•  • 

•■.00  At 

0.1-0  'A 

VI 

0,5'-V  x>.  •  . 

: 

.  1.  »  -■  * 

*'  -X'Ol 

O.  -  .f  V 

K’ -••••■O*.  0. 

.  ,  • 

.  • 

*  .0-  M 

.  «v  r.,i 

.  ,  A 

O .  *• 

A  .  0 A‘»  • 

'.A/,.1 ; 

0,r-*4  • 
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TABLE  A-7 


LISTING  OF  SECOND  API  PROGRAM 


VMOMIBPCD3V 

NOnISP 

f  i  :j  a  i  fem<-o 
F23  sr«-so 
133  or.  laifat+i 

I  ■’ll  t  M*  FM+  <  M*  +  /  <  -  Of  Nf  1-M  >  tS'L  )++/((  »M>-l>xMtST 
I  I  RTfSI+.xT 
I  A3  ->  <  A T  .'K  >  rGG 
I.  /  3  F  M 

v 

vHOMISRCtn* 

^  NOM1.SR 
FI  t  A  I 

r.'3  si«-sn 
r ‘  i  go :  at  fat+i 

I.  ‘13  EMf FM  LMx 1  -BTC 03 

cm  stfst+.xt 

C  A  3  ->(AT<K)pGO 

C  7  1  EM 

v 

vPERTARCim 
v  PERT  Aft 

0  3  '  PL  REFTT  TARGETING  RESULTS  FOR  ‘  S N *  *  TARGETS' 

>  l  1  8UBJEC1F  H  TO  1  >KF  SUCCESSIVE  ATT  ACK8 »  EACH  ATTACK  CONSISTING  ’ 

FO  1  OF  LIP  TO  •  JM»  '  MISSILES?  EACH  MISSILE  HAVING  A  PROBABILITY  OF  KILL  OF  ■  »Pf  '  .  1 

f  <13 

LSI  TM3 

f.  6  3  SCR-  (  NrO  )  »  1 

F?:.l  '  THE  EXPECTED  NUMBER  OF  MISSILES  FIRED  ISt  * 

[83  NOMISP 

C  9  3  'THE  INITIAL  STATE  ISt • 
r.103  SO 

rill  'THC  FINAL  STATE  IS S 1 
M?3  ST 

Cl.  33  'THE  SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  —  IT  ISt  ’  r+/ST 
M  *13  IH  IM+L  >-l 
'  IF. 3  EFT / 1 *ST 

'  1.6J  -I  HE.  EXPECTED  NUMBER  OF  TARGETS  SURVIVING  IS}  *  fE» 

I  1  VI  SD2M  +  /1  l*:»)xST)-<  I/T*ST)*2 
I  .18.3  Sli?FS»2rO 
!  :l  v  I  SDf-SD2#0  .  S 

C 7 0 3  '  THE  STANDARD  DEVIATION  OF  TARGETS  KILLED  ISt  ■ fSDJ 

v 

VRANTARCCI3V 
v  RANT  Aft 

II. 3  'RANDOM  I  ARRETING  RESULTS  FOR  '  ?Nr  '  TARGETS' 

123  '  SUBJECTED  TO  ' ? K 5 ‘  SUCCESSIVE  ATTACKSf  EACH  ATTACK  CONSISTING' 

r.Jj  '  OF  '  «Mf  MISSILES?  EACH  MISSILE  HAVING  A  PROBABILITY  OF  KILL  OF  'rP?'.' 

I  C.J  TM2 

II 6  3  SQM  NrO  >  1 1 

r  "7  I  *  T HP!  EXPECTED  NUMBER  OF  MISSILES  FIRED  ISt  ' 

VB3  NOMISR 

C  9  3  'THE  INITIAL  STAFF.'  ISt' 

11.03  SO 

Mil  'THE  FINAL  STATF  ISt' 

F 1 23  ST 

r>3J  'SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  —  IT  ISt  '  F  +  /ST 
I  M3  TFiN+1 
I  J“i3  Ff-F/i  x<5T 

Cl  6.1  '  THF  FXPFCTEP  NUMBER  OF  TARGETS  SURVIVING  IS?  *  ?E» 

L  1  J  3  SDX+/<I*2)xST>-<+/IxST>*2 

1 183  STi?PSB2rO 
C 193  SDfSD2*0.5 

C203  'THE  STANDARD  DEVIATION  OF  TARGETS  KILLED  IS»  *»SD» 
v 


A-14 


TABLE  A-7  (con't) 


VRANTARME039 
?  RANT  AW 

T 1 3  1 RANDOM  TARGETING  RESULTS  FOR  '»NS'  TARGETS  1 

L23  *  SUBJECTED  TO  *  fKi  •  SUCCESSIVE  ATTACKS?  EACH  ATTACK  CONSISTING ' 

€33  •  OF  • * H? •  MISSILES?  EACH  MISSILE.  HAVING  A  PROBABILITY  OF  KILL  OF  ' JP?  '  .  ' 

€43  1  (MISSILES  SPREAD  AMONG  CLOSEST  'fM»*  TARGETS)' 

€53 

CA3  TM1 
(  73  SCM-(NrO)  » 1 

€83  'THE  EXPECTED  NUMBER  OF  MISSILES  FIRED  ISt' 

€93  NOMISR 

003  'THE  INITIAL  STATE  IBS' 
fill  SO 

023  'THE  FINAL  STATE  IS?' 

033  ST 

0  43  *  SUM  OF  FINAL  STATE  ELEMENTS  SHOULD  BE  1  —  IT  XSS  '  H-/ST 

053  ION+1 
€.1A3  £«-  +  /!  xST 

073  'THE  EXPECTED  NUMBER  OE  TARGETS  SURVIVING  IS!  ' SE5 
L  3  83  SD2«-<+/<I*2)xST)-<+/IxST>*2 

093  SD?*SD2€0 

C203  SD«-SD2®0.5 

C21 3  'THE  STANDARD  DEVIATION  OF  TARGETS  KILLED  IS!  '  JSDJ 

■? 

vstiHcrm 
v  SIG«-SUM»D 
0  3  D«-~l 

C  23  SXG«-0 

€33  Cl<-J-I 
€43  EESD«-D+1 
C53  D1 8D-r  J 

C.A3  IGK-1*C1-D)X(DJC3  > x <Dl+< 1-Di ) XS)*M 
C73  SIG^-SIG+IG 
C83  -»(D<C1)pEE 

C93  STGOi;j)xSIG 

V 

VTM1C03V 
9  TfU 

0  3  0I0«-0 

T23  8M-P 

r:n  t«-<<n+i>?n+i>po 

i:  4 1  i«-~i 

E53  LOOPS  I<-X+l 

CA3  J«-“l 

r?3  lop:j«-j+i 

€83  -K((J-I)<0)v<J-I>>M)pBB 

€93  -KI=J)pCC 
0.03  -»(M<J)pYY 

013  ir  jsn^siGESUM 

023  BB S ■» ( J<N) pL.OP 
033  -»(I<N)pLOOP 
0  43  TC0»03<-1 

053  'THE  TRANSITION  MATRIX  T  IBS' 

063  '■ 

0  73  T 
003  ->0 

093  CC$T€I)J3«-S*M 

€20  X. — -9B® - 

C213  YY  S I  X«-I 
C223  JJ4-J 
€233  IM+M-J 
€24  3  J«-M 

€253  TC  JJf  I  I3*-SIG*-SUM 
C2A3  I«-II 
€273  J«-JJ 
C2B3  ->DB 


TABLE  A- 7  (con't) 


vihl'l  !. I  Iv 

V  TH? 

!  1  "!  fiiO'  0 

4  \  » 

C  3  3  1 t-  <  (  H  J-  J.  )  ?  r  H  t  !  •>  0 

r.  4  :i  i  «-“i 
C53  loop  :  y«-i+i 
1 6  :i  j*~j. 

171  LCif  t  J  1 

C  6  3  -><<<>- 1  >  <  0  >  v  <  j  - 1  )>M)p  BB 
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